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ABSTRACT 

In this paper we consider </> 2 scalar field potential as a candidate to dark matter. 
If it is an ultralight boson particle, it condensates like a Bose-Einstein system at very- 
early times and forms the basic structure of the Universe. Real scalar fields collapse in 
equilibrium configurations which oscillate in space-time (oscillatons).The cosmological 
behavior of the field equations are solved using the dynamical system formalism. We 
use the current cosmological parameters as constraints for the present value of the 
scalar field and reproduce the cosmological predictions of the standard ACDM model 
with this model. Therefore, scalar field dark matter seems to be a good alternative to 
cold dark matter nature. 
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1 INTRODUCTION 

Scalar fields are one of the most interesting and most myste- 
rious fields in theoretical physics. Fundamental scalar fields 
are needed in all unification's theories, however, there are 
not experimental evidence of its existence. From the stan- 
dard model of particles which needs the Higgs boson, until 
the superstring theory which contains the dilaton, passing 
throught the Kaluza-Klein and the Brans-Dicke theories or 
throught the inflationary model, scalar fields are necessary 
fields. Doubtless, if they exist, they have some features which 
make them very special. 

The Scalar Field Dark Matter (SFDM) model paradigm 
has been constructed step by step. One of the first sugges- 
tions that a (complex) scalar field could contri bute to struc- 
ture formation of the Universe was given bv IPressI (jl990h 
and lMadsenll| 19921 ). Nevertheless, complex scalar fiel ds were 
used b efore as matter candidates as boson stars bylftunini 
i|l969l ) (for a recent introduction to boson stars, see for ex- 
ample (GuzminF~s] (|2006h ). One of the first candidates to 
be scalar field dark matter is the axion, one of the solutions 
to the strong-CP problem in QCD (see an excellent review in 
iKolb. E. W. and Turner. S. T.I l|l99fj )'). Essentially, the ax- 
ion is a scalar field with mass restricted by observations to 
~ 10 -5 eV, which has its origin at 10 -30 seconds after the big 
bang, when the energy of the Universe was 10 12 GeV. This 
candidate is till now one of the most accepted candidates 
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for the nature of dark matter, if its abundance is about 10 9 
particles per cubic centimetre. 

The first in suggesting that a dark halo could be a Bose - 
Einstein condensate were [Sin| (|l994T ) and |ji fc Sinl (|l994T ) 
who used the weak field limit to show that a Bose-Einstein 
Condensate (BEC) with several nodes can fit the rotation 
galaxy curves with a very good accuracy . Further investiga - 
tions on this direction were performed bv lLee fc Kohl <|l996h . 
where they incorporated qb 4 interactions to the scalar field 
potential and used t he Gross- Pitaevskii equa tion instead of 
the Schrodinger one l|Leell 19961) . Nevertheless. ISeidel fc Suenl 
l|l99ll . I1994T ) showed that when the whole BEC is in the 
ground state, many nodes in Einstein-Klein-Gordon fields 
are unstable, since they evolve into the 0-node solution af- 
ter a while (for a clear explanation to this point see also 
i|Guzman. F. S. fc Urena-L6pezll2003l)) Thus, t he static so- 
lution s given by ISinl (|l994T l; |ji fc Sir] l| 19941 ); iLee fc Kohl 
(| 1996T ) are expected to be unstable. 

Later on, IPeebles fc Vilenkinl (|l999h proposed that a 
scalar field driven by inflation can behave as a perfect 
fluid and can have interesting observational consequences in 
structure formation. Besides that, they performed a sound 
waves analysis of this hypothesis giving some qualitative 
ideas for the evolution of th ese fields and called it fluid 
dark matter (IPeebl es 2000a, rj). I ndepe ndently and in an 
opposite way. iMatos fc Guzman! (1999) proposed a scalar 
field coming from some unify theory can condensate and 
collapse to form haloes of galaxies. Very early, this scalar 
field behaves as a perfect fluid, however its ultralight mass 
causes that the bosons condensate at very high tempera- 
ture and collapse in a very different way as the fluid dark 
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matter of iPeebles fc Vilenkinl ljl999t ) did. They were able 
to fit reasonably rotation curves of some galaxies using an 
exact solution of the Einstein equation s with an exponen- 
tial potential (Matos fc Guzman! 19991 ; iGuzman fc Matosl 
l200d ; lBernal. Matos fc Nunez ll200Sl ) . The first cosmological 
study o f the SFDM was performed in lMatos fc Urena-Lopea 
(2000313) 

where a cosh scalar field potential was used. The 
cosmology reproduces all features of the A Cold Dark Matter 
(ACDM) model in the linear regime of perturbations. 

On the other hand.lLesgourgues, Arbev fc Salatil (2002) 



and lArbev. Lesgourgues fc Salati ( 20031 ) used a complex 
scalar field with a quartic potential m 2 ^^ + X(cj>(j>^) 2 and 
solved perturbations equations (weak field limit approxima- 
tion) to fit the rotational curves of dwarf galaxies with a 
very good accuracy, provided that m 4 /A ~ 50 — 75 eV 4 . 

The importance of scalar fields in the dark sector 
has been increased, for instance, several authors have 
investigated the unification of dark matter and dark energy 
in a single scalar field ( Padmanabha n fc Choudhur 



2002; 



Arbev. Lesgourgues fc Salatil 



Bertacca. Barto lo & Matarrcsc 



cently 



Liddle fc Urena-Lopez 



2008). 



ihuryl 

200a ; 



Re- 
IHi); 

iLiddle. Cedric fc Urena-Lopez ( 200S ) proposed that the 
landscape of superstring theory can provide the Universe 
with a cf> 2 + A scalar field potential. Such scalar field can 
inflate the Universe during its early epoch, after that, the 
scalar field can decay into dark matter. The constant A can 
be interpreted as the cosmological one. This model could 
explain all unknown components of the Universe in a simple 
way. Another interesting model in order to explain the 
scalar fields unification, dark sector and inflation, is using 
a complex scalar field protected by an inter nal symmetry 
l|Perez-Lorenzana. Montesinos fc~M atos 200jl). 

In the present work the main idea is that if scalar fields 
are fundamental, they live as unified fields in some very early 
moment at the origin of the Universe. As the Universe ex- 
pands, the scalar fields cool together with the rest of the par- 
ticles until they decouple from the rest of the matter. After 
that, only the expansion of the Universe will keep cooling the 
scalar fields. If the scalar field fluctuation is under the criti- 
cal temperature of condensation, the object will collapse as 
a BEC. After inflation, primordial fluctuations cause that 
the scalar fields collapse and form haloes of galaxies and 
galaxy clusters. The cooling of scalar fields continue till the 
fluctuation separates from the expansion of the Universe. 

In this work we study the most simple model of SFDM, 
using a 4> 2 scalar field potential. In sections [2] and [3] we re- 
view the statistic of a boson gas to condensate and form a 
BEC, focusing in the necessary features for the BEC to form 
a halo of a galaxy and integrate the Einstein equations with 
a BEC matter. In section [4] we transform the Einstein field 
equations into a dynamical system, then we numerically in- 
tegrate them and look for the atractor points. We give some 
conditions on how these field equations can give the right 
behavior to reproduce the Universe we observed. Finally, in 
section [5] we conclude that this SFDM model could explain 
the dark matter of the Universe. 



2 THE STATISTIC OF A BEC 

In this section we review the condensation of an ideal Bose 
gas of N particles with mass m contained in a volume V 
with temperature T and with only a portion po of the sys- 
tem in the ground state. In order to see that and to be self 
contained, let us start from its grand partition function Q, 
which is given by 



Q(z,V,T) = Y[ Y~ze-^v' 



(1) 



where the fugacity z = e 13 ^ is defined in terms of the chem- 
ical potential p and /3 = 1/T. In this paper we use the 
fundamental constants % = c = kn = 1. 

Then, the equation of state for an ideal Bose gas is 

:logQ = -£>g(l-ze-^). (2) 
p 



T 



Thus, the grand partition function directly gives the pres- 
sure P as a function of z, V, and T. 

On the other hand the particle number N and the internal 
energy U are 

N = ^logQ = V , (3) 

az 1 — ze. pe p 



U 



d_ 

'df3 



-/3e p ■ 



(4) 



where e p is the single-particle energy with momentum p and 
the average occupation numbers < n p > are given by 



< rip >= 



1 - ze-^v ' 



(5) 



(6) 
(7) 



which satisfy the conditions 

N = J2 < n P >> 
p 

U = Yj e p < n p > . 
p 

Now we let V — + taking the limit of continuity, and 
replace sums over p by integrals over p, then we obtain the 
following equation of state 

py _ 2V r 

T - "(2^F7 

-/3p 2 /2m 



dp p 2 log(l — ze liv / 2m ) _ log(l — z), 



N 



2V /- 2 ze 
dp p 



(2^) 2 Jo 



1 - ze-PP 2 / 2 " 



+ 



1-z 



(8) 



These equations can be written into the equivalent form 



PV v r\ wi ^ 

— = ^3^5/20) -log(l-z), 



N 



V 



^3<?3/ 2 (X) + — 



(9) 
(10) 



where A = \J2-k jr-ciT is the thermal wavelength, and 

35/2(2) = ^= [ dx T 2 log(l - ze' 13 * 2 ), 

Jo 

53/2(2) = z—g 5/2 (z). 



(11) 
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Moreover, the internal energy is found from the formu- 
las ([2]) and © 



rr 3TV r ^ 

and as consequence the relation U=3/2PV is fulfilled. 
From equation ([5} we see that 



< n > = 



(12) 



(13) 



which is the average occupation number for a single parti- 
cle with occupation level p = 0. Equation (|10p can also be 
written as 



3 <no> 3 N , \ 

A — — = A — -gz/2{z). 



(14) 



This equation tell us that v > and therefore the 
temperature and the specific volume are such that X 3 y > 
53/2 i z )- This means that a finite fraction of the particles 
will be in the ground state with p = 0, i.e., the Bose 
gas condensates. In the region of condensation, the fugacity 
z ~ 1 and the function g(z) goes to the Riemann £ function 

sj(*W(0- 

The thermodynamical surface which separates the con- 
densation region from the rest of the P — V — T space, is 
given by 



a; 



N 
V 



= fls/a(l)= 2.612, 



(15) 



thus A c can be interpreted as the value for which the ther- 
mal wavelength is of the same order of magnitude as the 
average interparticle separation. Equation (I15|l defines the 
critical temperature for which the Bose Condensate forms. 
This temperature is given by 



T r _ = 



2tt 

5/3 



33/2(1) 



(16) 



where p — m^N /V is the density of the Bose gas. At con- 
stant temperature, equation (I16|l defines a critical density 



Pc 



«Vfl3/2(«) 



A :i 



(17) 



Thus, the region of condensation of the Boson gas is deter- 
mined by T < T c or p > p c . 

After the Bose gas condensates most of the bosons lie in 
the ground state, the scalar field starts to oscillate around 
the minimal of its potential and the scalar field starts to 
behave as dust (|Turner|[l983l) . Thus, after the scalar field 
decouples from the rest of the matter, the temperature of 
the BEC goes like 



Tbec = Tj 



(o) 



(18) 



where Tg EC ^ s * ne actual temperature of the BEC, a is the 
scale factor of the Universe and ao = 1 is the value of the 
scale factor at present. 

In the same way, as the BEC behaves as matter, its density 



goes like Pbec = Pbec/ 1 j where p BEC is the actual matter 
content of BEC in the Universe. With this result, equation 
(1161) can also be transformed into 



2tt 

573 



= 6.2 x 10 



where Qg EC is the actual rate of BEC, p cr it is the critical 
density of the Universe, h = _Ho/(100kms -1 Mpc -1 ) being 
Ho the actual value of Hubble's parameter. 

If the actual standard model of particles could be ex- 
tended to higher temperatures, we have to expect that the 
scalar field which forms the BEC, interacts with the rest of 
the particles to a temperature over some temperature T s . 
Because the physics of elemental particles is well known till 
temperatures like GeV, we do not expect that an exotic par- 
ticles as these scalar fields appear under temperatures like 
TeV. Here we have two possibilities, the first one is that 
the scalar field has never had interaction with the rest of 
the particles and it evolves independently from the rest of 
the fields, with only a gravitational interaction. In this case 
the scalar field condensates at the beginning of the Universe. 
The second possibility is that in the early universe the scalar 
field lived unified with the rest of the particles in a thermal 
bath and at some moment of its evolution, separates from 
the interaction. If this is the case let us suppose here that 
the scalar field which forms the BEC decouples from the rest 
of the matter at a temperature over TeV. Under this tem- 
perature, this scalar field has almost no interaction with the 
rest of the matter. If we expect that this scalar field forms a 
BEC, its critical temperature must be lower than the tem- 
perature of the scalar field decoupling. This fact gives us an 
upper bound of the mass of the scalar field 



O (0) n 
DEC i cri ^ 

C(3/2) 



- 31 (»i°Lft 2 ) 2/3 1 

(m^/GeV) 5 / 3 a 2 



GeV, 



(19) 



(20) 



m 4 < 10 -i 'eV. 



(21) 



On the other hand, from numerical simulations 
(|Seidel fc Suenlll99ll ) we know that scalar fields form gravi- 
tationally bounded objects with a critical mass given by 



(22) 



where m p i is the Planck mass and m is a factor such that 
m ~ 0.6 for both complex scalar fields (boson stars) and real 
scalar fields (oscillatons) . With the value given in (|2ip . the 
scalar field can form a gravitationally bounded BEC with a 
critical mass given by 



M crit > 1.491 x 10 64 GeV, 
= 2.658 x 10 40 gr, 
= 13.36 x 1O 6 M . 



(23) 
(24) 
(25) 



This is an interesting result, if there exists a scalar field 
and plays any role in the Universe at this moment, this scalar 
field must have a mass lower than the mass given in (|2ip and 
they are forming gravitati onally bounded BECs w ith masses 
around the mass given in lAlcubierre et alj (|2002l ). 



3 SELF-GRAVITATING BEC 

In this section we give some general features of the gravita- 
tional collapse of the BEC, we only pretend to show a generic 
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behavior of any self-gra vitating BEC. The BEC cosmol- 
ogy ha ve been studied bv lFukuvama. Masahiro fc Tatekawal 
l|2007l ) and many numeric a l simu l ations of this collapse are 



Alcubierre et al. d2002h; Guzman fc Urefia-Lopezl 
20041 ) ; [Guzman F. S. fc Urena-Lopez. L.I |2006l ) and be- 



sides. iGuzman. F. S. fc Urefia-Lopezl ( 2003 ) found that a 
BEC in the ground state are very stable under different ini- 
tial conditions. After the Bose gas condensates the gravita- 
tional force makes the gas collapse and form self-gravitating 
objects. Let us suppose that the halo is spherically symmet- 
ric, which could not be to far from the reality. In that case, 
the space-time metric reads 



ds 



2v 7.2 i 

-e at + 



dr 2 



+ r 2 dQ, 2 , 



(26) 



where the function v = v(r) is essentially the Newtonian 
potential and M = M(r) is the mass function given by 



Al 

dv 
dr 



= G 



4tt J p r dr, 
M + 4tt7- 3 P 



2MG \ 



(27) 



The Einstein field equations reduce to equations (|27 
and the Oppenheimer-Volkov equation 

dP _ (P + p)(M + 4-irr 3 P) 
~dr~ ~ r 2 (l- 2M£) ' 



(28) 



Let us focus in the case when the gas is far from forming 
a black hole. In that case we suppose that 2MG « r and 
equation (|28p reduces to 



dP 

dr 



= -4%Gr P{P + p). 



(29) 



The equation of state can be obtained from the equation 
PV = 2/3 U, JTOj and {T2). Combining all equations we 
obtain that 



2tt 



i>)2 j 53/2(2) 



gS/2(z) , .5/3 



= U>{p~ PO fl\ 

where u) is the constant 



2tt 95/2(2) 
m 8 / 3 ff3/2( ^ ) 5 / 3, 



(30) 
(31) 

(32) 



and po = rritt, < no > /V is the mean density of the par- 
ticles in the ground state. Thus, the Oppenheimer-Volkov 
equation (|28|) transforms into 



(33) 



g = -^Gr(p - p )(u(p - p ) 5/s + p). 



This differential equation can be easily numerically 
solved. Nevertheless, we have two interesting limits of equa- 
tion (|33p . First suppose that the uj constant is small such 
that P « p. This situation occurs for big scalar field masses 

~ rripianck- In that case, the equation (|33|l contains an 
analytical solution given by 

Po 



p(r) = 



I _ [ 1 £fi_ 

1 1 P(«) 



(34) 



where p(0) is the central density of the BEC. Observe that 




Figure 1. Plot of the p(r) function given in equation l|34j l for 
e < 1 (top plot) and for e > 1 (down plot). The plot is done in 
terms of p(r)/po- We have set e = 2 and e = 1/2 for each plot, 
respectively and po = 0.002. 



when r — ^ 00 , the function p(r) — 1 po- For numerical con- 
venience we set p(0) = epo in the plot, being e a constant. 
The function changes dramatically for different values of e. 
If £ > 1, the density p(r) decreases, but if e < 1 the density 
increases. The behavior of the density is shown in Fig. [1] 
This means that if the central density of the BEC is big- 
ger than the density of the ground state, we have the upper 
profile in Fig. [T] but if it is less than it, we have the bottom 
profile. 

The second and for us, a more interesting limit of equa- 
tion (|33|) is when P » p. This occurs when the scalar field 
mass is small enough m^, << mpianck, as for astrophysical 
BEC. In this limit the Oppenheimer-Volkov equation has 
also an analytical solution given by 



p(r) = 



P(0) - Po 



(2nGr 2 u>(p(0) - 
( P(0)/" 



po )5/3 + 1)3/5 

3/5 

+ P0, 



+ P0, 



(35) 



\2nGr 2 p(0) 

or equivalently P = l/(2nGr 2 + 1/P(0)). In this case the 
pressure dominates the BEC, the pressure acquire a maxi- 
mum for P(0)- Far away enough from the center of the BEC 
we can approximate equation (I35|l with 



(36) 



( 1/w \ 3/5 



which implies a space-time metric for the BEC given by 
dr 2 



ds 2 



1 - 2(r r 4 / 5 + UGp r 2 ) 



-exp(2v)dt +r 2 d£l z , (37) 
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the objects remain unalte red, since they behave like solitons 
l|Bernal fc Guzman! [2006). This means that in a merging of 
BEC they pass through each other without some alterations 
in its total mass as collisionless dark matter. This paradigm 
implies then that we must be able to see well formed galax- 
ies with the actual masses for very large redshifts, longer 
than those predicted by the bott om-up paradigm, i.e ., by 
CDM. In this sense some authors (|Cimatti et al.|[2004 ) sug- 
gest a discrepancy between the observed population of mas- 
sive spheroidal galaxies at high redshift with the numerical 
simulations of hierarchical merging in a ACDM scenario that 
underpredict this population. However, the discussion con- 
tinues because other physical processes, as feedback, could 
have important effects in this galaxies. 



10 



Figure 2. Rotation curve derived from metric (137 6 - The velocity 
and the coordinate r are in arbitrary units. 



where ro = 10/9(47r 2 /tj 3 ) 1/ ' 5 . Function v determines the 
circular velocity (the rotation curves) V ro t of test par- 
ticles around the BEC. Using the geodesic equation of 
metric (|37t one obtains that V 2 nt — rg t t,rf '(2<?tt) = ru' 
( Mat os. Guzman fc Nunez! I2OO0I ). Using equations (|27[l we 
can integrate the function v and obtain the rotation curves. 
The plot is shown in Fig. O where we see that the form 
of the rotation curves are analogous as the expected from 
t he observed in galaxies, specially in LSB and dwarf ones 
ilde Blok fc Bosmall2002l ;lde" Blok. Bosma fc McGaughll2003l ; 



ISimon et al.l 2*005h besides SFDM predicts a core den- 
sity profile that could have some astrophysics adva ntages 
l|Sanchez-Salcedo. Reves-Iturbide fe Hernandez! !200rJ l over 
the standard model (cuspy profiles). However, the discus- 
sion of the central region of the rotation curves continue. 
This is the main reason why it is not convenient to try self- 
gravitating BECs in the Newtonian limit. Remain that the 
Newton theory can be derived from the Einstein one for slow 
velocities, weak fields and pressures much smaller than the 
densities. However these last conditions is not fulfilled in 
self-gravitating BEC. 

From these results and from the simulations given 
m Guzman, F. S. fc Urena-Lopezl (|2003l ) it follows a novel 
paradigm for structure formation, which is different from 
the bottom-up one. In the SFDM paradigm, after the big 
bang the scalar field expands till decouples from the rest 
of the matter. If the scalar field has sufficient small mass 
such that its critical temperature of condensation is less 
than the temperature of decoupling, the scalar field forms a 
BEC. Then the scalar field collapses forming objects which 
final mass is not bigger than the critical mass m%i anck /m^,. 
These objects contain a density profile very similar to the 
profile shown in the top of Fig. [T] They are very stable un- 
der perturbations. It has been proposed that the dark mat- 
ter in gal axies and clusters is a scalar field with a mass of 
10" 22 eV (|Alcubierre et al.1 120021 ). If this were the case, the 
main difference for the structure formation of this ultralight 
scalar field with the bottom-up paradigm is that the SFDM 
objects form just after the collapse of the scalar field and 
remain so during the rest of the Universe expansion. Fur- 
thermore, they can collide together but after the collision 



4 THE COSMOLOGY 

In this section we review the Cosmology given by a 
SFDM model with two different scalar field potentials: 
V((f>) = \ rn 2 (j) 2 and V((f)) = Vb [cos!i(kA<^) — 1] where m is 
the mass of the boson particle, Vq and A are free parameters 
fixed with cosmological data and k, = 8tyG . Based on the 
curre nt observations of 5-year WMAP data (Hinshaw et al. 
2008) we will consider a Universe evolving in a spatially-flat 
Friedmann Lemartre-Robertson-Walker spacetime. We 
assume that this Universe contains a real scalar field (<j>) as 
dark matter, radiation (r), neutrinos (v), baryons (b) and a 
cosmological constant (A) as dark energy. 
The total energy density of a homogeneous scalar field is 
given by 



= 2p + v(<f>), 



the radiation and baryonic components are represented by 
perfect fluids with baryotropic equation of state p 7 = (7 — 
l)p 7 , where 7 is a constant, ^ 7 ^ 2. For example, for 



radiation and neutrinos (7r,i/ 



), for baryons (76 = 1) 



and finally for a cosmological constant (7A = 0). 

Thus, the field equations for a Universe with these com- 
ponents are given by 

2 

H = — ^{i 2 + iPi), 

4> + 3Hj> + d 4 ,V = 0, 
p 7 + 3-y H p 1 — 0, 

and the Friedmann equation 



(38) 



(39) 



In order to analyze the behavior of the different com- 
ponents of this Universe, we are going to use the dynamical 
system formalism following Appendix |A1 
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4.1 The <f> 2 scalar potential 



one. The whole close system is 



We start our cosmological analysis of SFDM taking the po- 
tential 



V(<t>)= 



(40) 



and developing the standard procedure to transform it into 
a dynamical system. For doing so, the new variables ()A2|) 
for the system of equations ((38} read 



K (f> 



k mtp 

~ V3H- <«> 
Using the definitions given in (|41|l . the evolution equa- 



tions (|38[) for potential (|40|) transform into an autonomous 
system 



/ 

X 


= — 3x — sit+^Ilx, 


(46a) 


/ 

u 


= sx+-llu, 


(46b) 




3 


(46c) 


I' 


= fni, 


(46d) 


1 

s 


- fn, 


(46e) 



In order to acquire geometrical information that dy- 
namical system analysis provide (see AppendixfX]), we study 
the stability of (I46p . To do this, we define the vector x = 
(x,u, z-,,1, s) and consider a linear perturbation of the form 
x — > x c + 8x. The linearized system reduces to 8x' — M8x, 
where M is the Jacobian matrix of x' and it reads as 



H 



o m 3 
— 3x — -j-m + -ila;, 
ri 2 

m 3_ 

I I+ 2 n "' 
3 

- (n-7) 2 7 , 

|(ar 2 +7*?) = ~n, 



(42) 



where as in AppendixfX] prime denotes a derivative with re- 
spect to the e- folding number = ln(a). Again the choice of 
phase-space variables (|41|l transforms the Friedmann equa- 
tion into a constraint equation 



F = x 2 + u 



= 1. 



(43) 



Because we are considering an expanding Universe 
which implies that H > and from the variable definitions 
(|41[) . we can see that u, z 7 ^ 0. With these variables, the 
density parameters can be written as 



Qdm = x 2 + u 2 , 

n, = 4 
ci a = i 2 , 



(44) 



where we have added explicity a cosmological constant vari- 
able I = za- Moreover, with the physical constraint SI ^ 
1 and the Friedmann equation Qdm + fi 7 + CIa = 1 the 
variable space is bounded by 

< x 2 + u + z 2 + I 2 < 1. 

On the other hand, observe that the variable space (|42p 
is not a completely autonomous one because H is an external 
parameter. In order to close the system we define a new 
variable s given by 



m 
77' 



(45) 



which dynamical equation (|A3d|) is 



s' = -lis. 

2 



M 



n - 3 + 6x 2 


—s 


37X 2 





— u 


6x u + s 


|n 


37 u z 





X 


Gx z 





§n + 372 2 -7 








6x I 





37 I z 


2 11 





6x s 





37 s z 





2 11 



The equilibrium points x 
{x, u, Zy, I, s}, considering only 7 



of the phase 
4/3, are then 



space 



(i) {±1,0,0,0,0} Kinetic scalar domination 

(ii) {0, 0, 1, 0, 0} Radiation domination 

(iii) {0,0,0, 1, s} Cosmological constant domination 

(iv) {0, u, 0,1,0} Cosmological constant and Potential 
scalar domination 

Finally, the eigenvalues of the matrix M valued at the crit- 
ical points listed above read 

(i) {6,3,3,3,3-7} 

(ii ) {i|,5a3 H+7 )) 



(iii) {0,0, |(- 

(iv) {-3,0,0,0, 



V9~ 

-7} 



-4s 2 ),i(-3 + V9-4s 2 ),-7} 



With this new variable, system (|42[) is now an autonomous 



As we can see, the radiation domination epoch shows 
a saddle point, however, in order to reproduce the big bang 
nucleosynthesis process is necessary that this kind of matter 
would had dominated in the past of the Universe. In other 
words, the radiation points should have corresponded to a 
source point. The domination of dark matter in the past (a 
source point) and the cosmological constant in the future 
(an attractor point) are showed in the Fig. 

In the following, we integrate system (|46[) with the con- 
straint (|43|) , following the procedure shown in Appendix 
[Al In general this system is very difficult to integrate be- 
cause it is a non-linear four-dimensional differential sys- 
tem of equations. It is clear that system f|46f) is a complete 
system which can fulfill or not the constraint (|43|l . How- 
ever, as it was shown in Appendix [X] system (|46[) together 
with constraint (|43|) is completely integrable. For simplic- 
ity we will take all the perfect fluid components as the 
equation z'^ — 3/2(11 — 7) z 7 with the Friedmann equation 
x 2 +u 2 + z 2 = 1. 
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Figure 3. Evolution of the density parameters for the system of 
equations l|46| l. The plot shows the dark matter domination epoch 
at early times, a source point. The cosmological constant at the 
future of the Universe is an attractor point. 



Thus, we substitute 3/2II from equation (|46e|l into the rest 
of the equations. With this substitution equation (|46c|) in- 
tegrates in terms of s as 



n^aexp(-- 7 JV), 



(47) 



where f2 7 is an integration constant. We multiply (|46a|l by 
2x and (I46b|) by 2u and sum both equations. We obtain 



(x 2 + u 2 )' = -6x 2 + 2 ln(s)'(x 2 + u 2 ). (48) 



Now, we use constraint (|43[) and equation (|47|) into equation 
(1481) to obtain 



6x 2 = 2 ln(s)' - 37 s 2 ^ 0) exp(-3 7 7V). 



(49) 



We substitute (gSJ) and (f4T)> into (|46ep to obtain 
= 0. Therefore, s is not an independent variable and 
we cast it into the system as a control variable which 
parametrizes the decrease of H, a s imilar result is found by 
lUrena-Lopez fc Reves-Ibarral l|2007h . In what follows we will 
use this important result. 

Of course, to guess variable s in order to fulfill con- 
straint (143[) is not so easy. In order to avoid this problem we 
can consider the observed dynamic for H and model it by 
the following ansatz 



H 



to* 



(50) 



because it is well-know the behavior for H at different epochs 



2? ^ 



(51) 



There exists a restriction in the parameter n. Because is 
well know that H is a function monotonically decreasing, n 
has to satisfy n ^ 0. With the ansatz l|50[l . the dynamical 
equation for s reads 



s' = (mt )~ 



11 I - 

s 



so s 



(52) 



where we have defined k = 1/n — 2. 

In the following, we investigate if this system can repro- 
duce the observed Universe. We introduce the components 



of the background Universe into the dynamical system de- 
scribed by (|46p adding to it baryons (b), radiation (2) and 
neutrinos (v). Thus, the system transforms into 



X 


= — 3 x — su -\- — fix, 


(53a) 


u 


3 

= sx+-Tlu, 


(53b) 


b' 


= ^(n-i)6, 


(53c) 


z 


- IH)- 


(53d) 


V 


- §H), 


(53e) 


I' 


, |n, 


(53f) 


s 


= S S~ , 


(53g) 



with II = 2x 2 + b 2 + |z 2 + and the Friedmann equation 
reduces to the constraint 



F = x 2 + u + b 2 + z 2 + v 2 + I 2 = 1. 



(54) 



Using this ansatz we can reduce till quadratures the solution 
of system (|53[) . In order to do this, observe that 

^rr -fc-l 

-n = s s 

Now, using this last identity, equation (|53c[l - (|53f|l can be 
integrated to give 



z 1 = 20 [s (k + 1) N + si] *+i e" 



for each corresponding value of 7. Finally, equations (|53a|l 
and (|53b|) can be integrated as follows. We divide (|53a[l by x 
and (|53b|) by u and take the difference between both equa- 
tions. We define y = x/u to obtain 

y' + 3y + q(N)y 2 = -q(N), (55) 

where function q(N) = [so (k + 1) N + Equation 
(|55|) is a Riccati equation which can be reduce to a Bernoulli 
equation by defining y — w+yi, where yi is a known solution 
of (1551). It reduces to 



It/ + (3 + 2qy±) w + q z 2 = 0. 



(56) 



Equation (|56|) can be further reduced by defining W = 1/w, 
we obtain 



W' - (3 + 2gyi)IU-g = 0, 
which integral is 

W = e A J e~ A qdN, 
with A = / (3 + 2 q j/i) dN. Thus 
u = uq q exp I yq dN ) , 



-3JV 

x = xoqe exp 



q -dN 

V 



zoqe 2 
1 



7 JV 



W 



+ yi- 



(57) 
(58) 

(59a) 

(59b) 

(59c) 
(59d) 



In the particular case where so = 0, the integrals can be 
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Figure 4. Evolution of the density parameters for the system 
(1 53 I t with n = 1 (top panel) and n = 5 (bottom panel). This 
values of n are not reproduce the standard behavior of ACDM 



solved analytically, however this value for so does not have 
a physical meaning. 

On the other hand, we can evaluate the integrals us- 
ing numerical methods for different values of the free con- 
stants. We can obtain a numerical solution for the sys- 
tem using (|59|l or directly integrating system (|53|l with an 
Adams-Bashforth-Moulton (ABM) method and using as ini- 
tial data the WMAP+BAO+SN recommended values to 
f^ 0) = 0.721, VL ( ° ] M = 0.233, n ( b 0) = 0.0454, n ( r 0) = 0.0004, 

= 0.0002, the result is the same. 

Fig. [4] and Fig. [5] show the numerical solutions of the 
dynamical system (I53|l . In Fig. 2] we set n 1, as exam- 
ples we show n = 1, 5. From these figures it is clear that 
the radiation remains subdominant for this values of n. On 
the other hand, in Fig. [S] where the plots were made for 
n — 1/2, 1/5, the radiation and the neutrinos behave ex- 
actly in the same way as they do in the ACDM model so we 
expect that both of these can reproduce the observed Uni- 
verse. The first values for n are not able to explain the big 
bang nucleosynthesis, since radiation never dominates as it 
is required. However, the last values for n can reproduce the 
radiation dominated era. Following the radiation dominated 
era, <j) 2 dark matter becomes the component that dominates 
the evolution and finally the Universe is dominated by the 
cosmological constant. Fig [6] shows the constraint F in (|54[) 
in order to visualize the integration's error. Observe that 
F w 1 at every point in the evolution, indicating that the 
Friedmann equation is exactly fulfilled all the time, this be- 
havior is exactly the same for all runs. 
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Figure 5. Upper panel: evolution of the density parameters for 
the system (1 53 I t with n = 1/2. Lower panel: evolution of the 
density parameters for the system H53H with n = 1/5. SFDM 
reproduces the standard ACDM behavior in both cases 
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Figure 6. Evolution of the function F = x 2 +u 2 +b 2 + z 2 + v 2 +l 2 
in H54H for the system J53H with n = 1, 5, 1/2 and 1/5. Function 
F is exactly the same for all values of n in all these cases. 



4.2 The cosh scalar potential 

Now, we are going to compare above results with the poten- 
tial 

V{<f>) = V [cosh( K A<^) - 1] . (60) 
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In order to do so, we define new variables as 

= 7Eh> 



2 Vq k 
3 



cosh 



2V k 
3 H" 



sinh 



I ; 



k\4> 



VSH^-VSH- (61) 
Substituting definitions (|6l) into equations ((38|) we ob- 



tain 



—3 x — \vu + -III, 

AXV + —II M, 

3^ 
Xxu + — II w, 



i' = 



(n- 7 ) »y, 



(62) 



where again the prime means derivatives with respect the 
e-folding number N = ln(a) and we are also using the func- 
tion IT = 2a; 2 + 7Z 2 . From the definitions (|61|l it follows the 
constraints 

i 2_ )2 = 2V k 2 1 _ J_m| 
A 2 # 2 ' 



(63) 



3 H 2 

and the Friedmann equation (|43[) written in this variables 
reads 



2 , 2 . 2 i t2 
X + M + Z + I 



(64) 

However, equation (I64p is actually not real constraint, 
since they are inhered in the dynamical equations (|62[1 (see 
appendix A equation (IA6|l '). Furthermore, constraint (|63p 
is also inhered in the dynamical system, observe that if wc 
multiply the second equation of (|62|l by 1/2 u and the third 
by 1/2 v and rest each other, we obtain 



(65) 



But this relation follows directly from the field equa- 
tions (|38[) . This means that system (162[) is compatible with 
the constraint (|fc>3[) . Using this constraint (I63|) in the dynam- 
ical system (I62|l . we obtain 



-3 x - u^A^+^j 1 + : 



/ ra \ 2 3 
= :r^/A 2 u 2 + ^— J + -Hu, 



(66) 



/' = in, 

We notice, that occurs the same situation as (f> 2 po- 
tential. Introducing again the variable s = m/H with its 
dynamical equation. 



s' — (mto) n 1 n ( — 



(67) 




Figure 7. Evolution of the density parameters for the system 
l|68j l. where the scalar field potential is given by the equation l|60|l . 



we obtain 



-3 x — u\J \ 2 u 2 + s 2 + — Ila;, 



I' = 



xsj X 2 u 2 + s 2 + -IIu, 

1 (n- 7 ) a,, 

-Hi, 

2 ' 



•so 



(68) 



The density parameters are the same as we have defined 
at (|44[) . We solve numerically (|68|l with the same initial con- 
ditions as the system of equations (|53p and A w 20. The so- 
lutions are shown in Fig. (0 ■ The plot shows the dynamical 
evolution for a Universe with SFDM with the potential (|6Up . 
notice that is equivalent to potential (|40[) . 

Finally, we use the same dynamical system formalism 
for the case of ACDM in order to compare with SFDM. 
We consider that it background Universe is composed by 
baryons, radiation, neutrinos, cold dark matter and cosmo- 
logical constant with an equation of state as perfect fluid. 
We solve numerically this system and in general terms the 
dynamic of both scalar potentials is indistinguishable of the 
standard model. This is an important goal of this paper. 

The next step is to compute the age of the Universe 
using our model. The age equation can be written as 



to 



(69) 



Using the definition for I from (|4T1) or ((61) . eq. (f69l) 

reduces to 



to = 



K V P A J No 



IdN. 



(70) 



We compute (JTOJ) and obtain that t a ~ 13.77 Gyr. This 
result is in agreement with the cosmological observations 
from WMAP+BAO+SN which estimate t„ = 13.73 ± 0.12 
Gyr and therefore H a = 70.1 ± 1.3 kms~ 1 Mpc~ 1 . Further- 
more, in Fig. [7] we see that scale factor of decoupling is 

10" 3 , this means a redshift z ~ 1000. At this redshift, 
the neutrinos made up ~ 12% of the Universe. On the other 
hand, WMAP cosmological observations show that when the 
Universe was only 380,000 years old, neutrinos permeate the 
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Universe within 10% of its total energy density. Thus, SFDM 
is in agreement within the measurements of WMAP. This 
result shows that scalar field is a plausible candidate for dark 
matter because it behaves like cold dark matter. 



5 CONCLUSIONS 

SFDM has provided to be an alternative model for the dark 
matter nature of the Universe. We have shown that the 
scalar field with a ultralight mass condensates very early 
in the Universe and generically form BEC's with a density 
profile which is very similar as that of the CDM model, but 
with a almost fiat central density profile, as it seems to be in 
LSB and dwarf galaxies. This fact can be a crucial difference 
between both models. If the flat central density is no con- 
firmed in galaxies, we can rule out the SFDM model, but if 
this observation is confirmed, this can be a point in favor of 
the SFDM model. We also show that the l/2m 2 ^ 2 potential 
and the Vo[cosh(«;A<jf>) — 1] model are in fact the same. They 
have the same predictions, a control variable which deter- 
mines the behavior of the model, given naturally the right 
expected cosmology and the same cosmology as the CDM 
model. This implies that the differences between both mod- 
els, the CDM and SFDM ones, is in the non linear regime 
of perturbations. In this way they form galaxies and galaxy 
clusters, specially in the center of galaxies where the SFDM 
model predicts a flat density profile. If the existence of su- 
persymmetry is confirmed, the DM supersymmetric parti- 
cles would be observed by detectors and they would have 
the right mass, DM density and coupling constant, there- 
fore the SFDM model can be ruled out. However, if these 
observations are not confirmed, the SFDM is an excellent 
alternative candidate to be the nature of the DM of the 
Universe. 
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x e X, 



(Al) 



where f : X — * X . 

The main step to get qualitative information on 
solutions is studying the flow of the equation in the vicinity 
of their critical points based on the Hartman-Grobman 
theorem, namely the study of its stability. 

The essential idea is firstly find the fixed (or critical) 
points of the equation (|A1|) which are given by f(x c ) = 0. 
Then linearized the differential equation at each critical 
point, that is, expanding about these points x = x c + 5x 
which yields to 

5x = MS5c, 

where M is the Jacobian matrix of x'. Therefore the general 
solutions for the linear perturbation evolution can be written 
as 

r— / r~ Af<5t 

ox = oxoe , 

where N is the matrix composed of the eignvalues rm asso- 
ciated to AL 

The stability of the system (|A1|I depends on the values 
of the eigenvalues: if the real part of all eigenvalues is nega- 
tive, the fixed point is asymptotically stable, i.e., an attrac- 
tor. All eigenvalues with positive real part make the fixed 
point asymptotically unstable (commonly called as source 
or repeller). 

On the other hand, a saddle point happens when there 
exists a combination s of stable a nd unstable points. For a 
extended review see, IColevi i| 19991 ). 

Then, we give a procedure for transforming equations ((38} 
and (|39[) . with an arbitrary potential, into a dynamical sys- 
tem. We define the dimensionless variables 



K (f) k VV 

x = V6~h> u = 73"r' 



V3 H 



(A2) 



Using above definitions (|A2|I . the evolution equations 
(1381) transform into an autonomous system 



H 



-(2x 2 +7Z 7 ) 



-If. 



(A3a) 
(A3b) 
(A3c) 
(A3d) 



APPENDIX A: 
REVIEW. 



DYNAMICAL SYSTEM 



The theory of dynamical systems is used in the study of 
physical systems that evolve over time. It is assumed that 
the physical state of the system to an instant of time t is 
described by an element a; of a space phase X, which can be 
of finite or infinite dimension. The evolution of the system 
is represented by a differential autonomous equation in X, 
written symbolically as 



This last equation (|A3d[) can be written also as 



a =-n s , 



(A4) 



for the variable s — cte./H , and determines the evolution of 
the horizon. Here a prime denotes a derivative with respect 
to the e- folding number N = ln(a) . The Friedmann equation 
(|43|) transforms into a constraint equation 



2 . 2.2 

X + u + z~ 



1. 



(A5) 
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With these variables, the SFDM density can be written as 

Om/ = x 2 + u 2 . 

Observe that if we derive (|A5|l with respect to JV and 
substitute system (|A3f) into this, we obtain 

F' = 3(F-1)II, (A6) 

indicating that constraint (|A5|) is compatible with system 
(|A3|) for all scalar field potentials if the Friedmann equation 
is fulfilled. 

Now we show that system <|A3[) together with constraint 
(IA5I) is completely integrable. To integrate system (|A3[l . first 
observe that we can substitute 3/211 from equation (|A4|) into 
the rest of the equations. With this substitution equation 
(|A3c[) can be integrated in terms of s as 

z 7 = 1/0^ sexp(-^N), (A7) 

where 0.^ is an integrations constant. Now we multiply 
l!A3al) by 2 a; and (|A3b|) by 2 u and sum both equations. We 
obtain 

(x 2 + u 2 )' = -6x 2 + 2 ln(s)'(:r 2 + u 2 ). (A8) 

Now, we use constraint ()A5[) and equation (|A7I) into 
equation (|A8|) to obtain 

6 a; 2 = 2 ln(s)' -37s 2 ^ 0) exp(-3 7A^). (A9) 

Now we have to integrate equation (|A4I) with all these 
results. If we substitute (|A"9)| and (|A"7)) into (lA3dl) or (|A"4| 
we obtain = 0, that means s is an arbitrary variable which 
parametrizes the decrease of H and can be cast into the 
system as a control variable. In other words, equations (|A3d|) 
and (|A4|) are actually identities, and not equations. 

Thus, we set the variable s from system (I A3f) arbitrary 
in the equations (|A3a|l . (IA3b|) and (|A3c[) . 
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